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I We investigate antiferromagnetic spin ladders with nonmagnetic impuri- 

O 

\^ ' ties by variational and numerical (Lanczos and DMRG) methods. The inter- 

P\ ' 

^ action between the two unpaired spins opposite to the impurities is described 

S I by an effective exchange interaction Jeff, the magnitude of which depends 

^ ' on the impurity distance. The magnitude of Jeff is different for unpaired 

o ■ 

J-^ ' spins at the edges of an open ladder and in the bulk. This difference is le- 

lated to the different distribution of the unpaired spin into the bulk of the 

' ladder. The numerical results are interpreted using matrix product states. 

Using the DMRG we calculate the spectrum of low-lying energy levels for up 

to 6 impurities and find that these spectra can be reproduced assuming pair 

interactions with an accuracy of better than 10%. We discuss the filling of the 

ladder gap with impurity states and argue that in the thermodynamic limit 

the spin ladder with a finite concentration of impurities always shows a Curie 

susceptibility at low temperatures. 
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I. INTRODUCTION 



In the last few years spin systems consisting of two interacting spin chains with S = 
|, now usually called (two-legged) spin ladders have attracted considerable attention as 
recently reviewed by Dagotto and Rice§. The spin ladder with isotropic antiferromagnetic 
interactions exhibits a spin gap in the excitation spectrum which has been shown to be 
related to the dimer gap (af coupling on the rungs only) as well as to the Haldane gap 
(strong ferromagnetic coupling on the rungs). It has been concluded! that the isotropic spin 
ladder is in the same phase as the Haldane chain - this can formally be described using 
matrix product ground states in a generalized spin ladder with additional diagonal bond^. 

In the Haldane chain, the investigation of magnetic as well as nonmagnetic impurities has 
contributed substantially to our understanding of the system: There exist quasifree spins 
5 = ^ at the end of the chain segments created by the impurities as qualitatively predicted 
by the valence bond ground stateB and experimentally found in ESR experiments on NENP 
doped with cJ. We have investigated nonmagnetic impurities in spin ladders and found 
that these lead to effects both similar to and different from impurities in the Haldane chain. 
The essential difference is of geometric origin: impurities in a spin ladder do not break the 
sequence of magnetic interactions, but create weak links: one spin S = ^ on one leg of 
the ladder has no counterpart on the other leg and mediates between the adjacent regular 
ladder structures. In the following we use the term 'unpaired spin' for these spins without 
counterpart on the opposite leg. A typical configuration with three impurities is shown in 
fig. |l|; neighboring impurities may lead to unpaired spins on the same leg (cis configuration) 
or on opposite legs (trans configuration). This structural effect that the impurities do not 
break the coupling along the linear arrangement is analogous to the one of impurities in 
a chain with nn and nnn interactions - which can actually be considered as a generalized 
ladder!. For a statistical distribution of impurities both the sequence of cis and trans pairs 
and the distances between these pairs will be statistical. The distance between neighboring 
unpaired spins is conveniently counted by the number p of complete rungs between them. 
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Assuming small impurity concentration c we neglect the possibility of complete breaking of 
the ladder which is proportional (?. 

For a microscopic understanding of the effects of a low concentration of defects on the 
spectrum of the spin ladder we present in the following numerical and analytical calculations 
for the properties of the ground state and of low-lying excited states of impure ladders with 
5 = 1 and all exchange interactions antiferromagnetic, isotropic and of equal magnitude 
(which is set equal to unity). In particular we have considered: 

• two spins 5* = I in either cis or trans configuration at the edges of an open ladder 
formed by p rungs (exact diagonalization with the Lanczos algorithm and analytical 
estimates), 

• two impurities at varying distances in either cis or trans configuration in the bulk ladder 
(exact diagonalization with the Lanzcos algorithm for periodic boundary conditions, 
DMRG calculations for open boundary conditions and estimates from matrix product 
states), 

• the spin ladder with many impurities (up to six impurities in DMRG calculations and 
arguments based on the Lieb-Mattis theorem). 

Our work is motivated by recent experimental on the ladder material SrCu203 doped 
with Zn, which show a tendency of the ladder gap to vanish with increasing Zn concentration. 
Calculations for similar systems which, however, concentrate on somewhat different aspects 
have recently been done by Motome et al0 and Martins et all. 

II. NUMERICAL METHODS 

We have performed exact diagonalizations using the Lanczos technique for spin ladder 
systems with two defects and a varying number p < 12 of complete rungs, i.e. N = 2p + 2 
spins in total and have considered both periodic and open boundary conditions. For open 
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boundary conditions we study a ladder with two impurities at the edges of the system and 
consequently separated by p rungs in either cis or trans configuration. In the case of periodic 
boundary conditions we fix one defect on the first rung and vary the distance (pi, pi < [p/2] 
rungs) to the second defect. We again consider cis and trans configurations with p = 11 and 
p = 12 rungs. For these configurations we have calculated the energies of the ground state 
and of low lying excitations and, for open boundary conditions only, the distribution of the 
magnetization. 

For the normal ladder system projection of the Hilbert space on the irreducible rep- 
resentation of the corresponding symmetry groups (refiection and mirror symmetry, total 
SU(2) invariance and, in the case of periodic boundary conditions, translational symmetry) 
drastically reduces the amount of memory for diagonalization within these subspaces. As 
a conseqence of inserting impurities into the ladder most of these symmetries, in particular 
translational invariance, are no longer present and the computational effort is correspond- 
ingly larger. It is only the SU(2) invariance of the model which is preserved. 

The calculations were performed on a MPP CRAY T3D SC256 of the Zuse Computing 
Centre Berlin using a parallel implementation of the Lanczos algorithm. The numerical 
accuracy of these calculations is 10^^° or better. 

Whereas the accuracy of the exact diagonalizations using the Lanzcos approach is high, 
it can only be done for relatively small systems. To study longer chains we have used the 
Density Matrix Renormalization Group (DMRG)i. Using this method we have studied open 
ladders from 2 x 50 up to 2 x 100 spins. As the defects are found to behave like spins localized 
on the bulk correlation length scale, defects more than 20 sites away from the ladder ends 
behave effectively like in an infinite system. 

As good quantum number we used the total spin; in configurations with defects 
distributed symmetrically around the ladder center, we also used parity. This allows for fast 
classification of states. 

We typically kept M = 100 to M = 150 block states. A remark is in order on the precision 
of the DMRG in this particular application. For the defect-free spin ladder, truncation errors 
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are very low, for M = 100 p = 2.9 x 10~^^ and for M = 150 p = 1.1 x 10~^^, corresponding 
to an extremely small error in the energies of the defect-free system for low-energy states. 
The effective precision of the DMRG is however greatly reduced, when a defect or defect 
pair are added. The block states from the step before are not optimally chosen to represent 
the altered system. We find that the introduction of defects effectively reduces the precision 
of the DMRG by up to several orders of magnitude. 

This problem becomes particularly pressing in defect configurations nonsymmetric with 
respect to the ladder center. For any even number of defects, the ground state is the lowest 
eigen state in the 5*^ = sector. A nonsymmetric defect configuration implies that during 
the growth process there will be an odd number of defects present in some DMRG steps. 
The total number of spins on the ladder changes from even to odd. Consequently, the 
ground state jumps to the = ±| sector, to which the available block states are not well 
adapted, drastically reducing the precision of the DMRG. However, the error can be greatly 
reduced by using the finite size DMRG algorithm^: The whole chain is recalculated in 
the presence of the complete defect configuration without any jumps between different 
sectors. The gain in precision by far exceeds the one obtained by increasing M. 

To illustrate this point, let us consider the case of a 2 x 50 spin ladder with two defects 
in trans-configuration on rungs 24 and 26. We have calculated the lowest eigenstates with 

= and 5*^ = 1, keeping M = 40 and M = 100 states. In Table | the energies given are 
for the chain calculated by conventional DMRG, and after one resp. two iterations applying 
the finite size algorithm. Obviously, even for M = 100 the unmodified DMRG (first line) 
gives results of the order of 0.01 away from the converged result, whereas the finite size 
algorithm produces highly precise results even for M = 40 (compare results in second and 
third lines). 

All DMRG calculations were performed on a PentiumPro 200MHz machine running 
under Linux. 
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III. TWO UNPAIRED SPINS 



In this section we discuss the spectra and spin configurations for two unpaired spins in 
ladders with both open and periodic boundary conditions, presenting and comparing nu- 
merical results from both Lanczos and DMRG calculations and from analytical approaches. 

Numerical results for the spectrum of unpaired edge spins 

We consider the configuration with p rungs connecting two spins S, S' in either cis or 
trans configuration at the ends. Using the Lanczos algorithm we have exactly diagonalized 
systems with p < 12. We find that there is an energy gap of the order of the known ladder 
gap (A ^ 0.53) but the ground state is splitted into a singlet and a triplet state. As a 
consequence of the Lieb-Mattis theorem the lowest state for our configuration is a singlet if 
the two unpaired spins are on different sublattices {p even for cis spins and p odd for trans 
spins) and a triplet if they are on the same sublattice {p odd for cis spins and p even for 
trans spins) as discussed in ref. |12|. We interpret this splitting as resulting from an effective 
interaction between the unpaired spins 5*, S' at the ladder boundaries and we derive an 
effective coupling by writing a Hamiltonian in the subspace of these lowest two states as 

H^I^ = Eo,^ +jf/,%ip)SS'. (1) 

The index ± refers to the sign of Jeff, i-e. antiferromagnetic and ferromagnetic effective 
interaction respectively. J^ff'^'^ip) is given in Table |T| and is also plotted in fig. |^ to show 
its dependence on p. Jeff% actually is identical to the singlet-triplet splitting of the ground 
state. An excellent fit for p > 5 is obtained as 

Eo,+ (p) = -1.157 

J:;f±(p) = Jo,± e-^/«, with Jo,+ ~ 0.674, Jq,- ~ 0.714, ^ ^ 3.1. (2) 

The uncertainty in these data as determined from a least square fit is 10~^ for Jo,± and better 
for -Eo,+- Due to the finite size of the open chain, i?o,- (~ —1.161 p) cannot be determined 
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with equal accuracy. We draw attention to the fact that the best fit is purely exponentially 
decaying and that a behavior oc exp{—p/^)/y/p can be excluded. This is in parallel to the 
behavior of an open Haldane chaini. 



Spectrum of unpaired spins in the bulk ladder: Numerical results 

Contrary to the Lanczos approach the DMRG allows to deal with sufficiently long ladders 
so that the two unpaired spins can be at some rather large distance and at the same time 
sufficiently far from the boundaries to identify their bulk interaction. The splitting of the 
ground state for two unpaired spins is again described by a Hamiltonian as above with 
Jeff\{p) replaced by Jgjj-\(p). The latter quantity is also given in Table || and shown in 
fig. We find (for p > 5) 

J^^ff%{p) = J^l^e~P/^, with Jo^^^f'^ ^ 0.43, e~3.1. (3) 

Jq^'^ is approximately equal to Jo,+^' ^ut determined less accurately. We notice that J^jfip) 

'Age 
ff 

Jeff%{p) =re.cAP)J:ff:±{p)- (4) 



is reduced with respect to J^fTip) by a reduction factor rexc,± 



For effectively antiferromagnetic interaction the reduction factor rg^c is p independent, 
i^exc ~ 0.65, for effectively ferromagnetic interaction it approaches this value asymptotically. 
This reduction will be explained quantitatively below; qualitatively it is due to the fact that 
the unpaired spin projection delocalizes both into its right and left neighborhood in the bulk 
case whereas there are neighbors on one side only in the boundary case. The dependence 
of the effective exchange on distance is again purely exponential and the prediction in ref. 



T2| , which includes an additional factor l/i/p is not verified from our data. The 'best' 
fit enforcing the factor l/^^P ends up with ,^ ^ 4, i.e. a correlation length which is not 
appropriate. 

Using the Lanczos algorithm, we have also done exact diagonalizations for configurations 
with periodic boundary conditions (pbc) and two spins S, S' in cis or trans configuration 



separated by pi, resp. p2 rungs. The spectra are shown in fig. for pi + p2 = 12 and 
in fig. ^,d for pi + P2 = H; their quahtative structure is the same as for the open ladder. 
We find that the smaller of the two values pi,p2 determines whether the ground state is 
singlet or triplet when the rules given above for the open ladder are used. From the graphs 
presented in fig. ^ it is also evident that a much larger singlet-triplet splitting of the ground 
state is obtained for pi + p2 even than for pi + p2 odd. The obvious reason is that the two 
exchange interactions which have to be added for pbc have equal (different) sign for pi +p2 
even (odd). Using again the concept of an effective exchange interaction, Jef'f{pi,P2), we 
expect for two unpaired spins (resulting from two impurities in a periodic ladder) 

jtf}{PuP2) = j'effiPl) + j'eff{P2). (5) 

Inserting the results as given before, we find that this relation is obeyed to within 10%. 

In fig. ^ we also show the results for the lowest excited states, i.e. in particular the 
splitting of the regular ladder gap with energy ^ A. For most of the defect configurations 
the effective interactions have different sign for the ground state and for the first excited 
state (energy of the regular ladder gap). For example, a singlet ground state is related to 
a triplet as the lower one of the states of energy A, as follows from a simple coupling of 
the impurities to the lowest bulk excitation with S = 1. This rule, however, is not strictly 
obeyed and a detailed discussion of the structure of the excited states in the impure ladder 
has to be reserved for future work. 

Analytical approaches to unpaired spins on the ladder edges 

A qualitative understanding of our numerical results for the open ladder can be obtained 
using matrix product wave functions as described in ref. |^. We consider a system with two 
spins at the ladder edges in trans configuration and p rungs in between; A wave function 
describing this system with four degrees of freedom can be written down as a 2 x 2 matrix 
product (MP) wave function. 
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p+1 



9j 



( a\tQ)j+h\s)j -aV2|t+), ^ 



(6) 



Here the matrix Qj describes the couphng of spins situated on diagonal sites on two adjacent 
rungs to singlets, resp. triplets; the boundary spins are coupled likewise to the neighboring 
rung. The value of 6^ = 1 — 3a^ is determined by minimization of the energy; for the 
numerical estimates below we take the result for the infinite laddeJi, h ~ 0.1735. We have 
chosen the two edge spins in trans configuration where a MP wave function is easily written 
down, whereas this is more difficult for two spins in cis configuration (with e.g. two more 
spins on the upper leg than on the lower leg); the latter case requires a more detailed analysis 
which will be published separately. Instructive limiting cases of the wavefunction of eq.@ 
are: 

• For a = 6 = I it is the wave function for p singlets on the p rungs and truly free bound- 
ary spins 5, 5". This wave function actually is an eigenfunction to the Majumdar- 
Ghosh Hamiltonian for a finite chain with 2p spins and one additional free spin at 
each open end. Their spin projections can be identified with the matrix indices cr, a' . 

• For b = 0,a = l/^/S, |V')o-,cr' describes p+1 units with 5 = 1 on the diagonals of 
the ladder, i.e. it is identical to the 4 eigenfunctions for the AKLT chain with open 
endsi. In this limit it is known that at the edges there exist quasifree S = ^ spins 
which extend somewhat into the bulk of the ladder. Note that the matrix indices a, a' 
now can no more be identified with the spin projections of the edge spins. 

• The general case will be intermediate and the boundary spins extend into the adjacent 
bulk ladder to an extent determined by the correlation length. This correlation length 
is underestimated by the MP ansatz but the tendency of the variation is given correctly: 
it decreases with increasing b, i.e. from the AKLT limit via the ladder configuration 
to the Majumdar Ghosh limit. 
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For a more quantitative treatment we couple the 4 MP wavefunctions Hi 9 ? cr? cr' = 
±1, to singlet, resp. triplet states and calculate the effective exchange energy from the 
energy difference of these states. This leads to the following results for a configuration with 
p complete rungs in between: 

jfff':Mp{p) = ±Joe-''/^, with Jo ^0.77, r'~l-23. (7) 

This gives the correct sign of Jg// (i.e. ferro- or antiferromagnetic if a, a' are on the same 
or on different sublattices respectively) and a remarkably good numerical value for p = 0; 
the decay with distance, however, is too strong, due to the fact that the correlation length 
is underestimated in MP states. 

We notice that the low energy structure of the spectrum is the generalization to spin 
ladders of the Kennedy triplet observed in Haldane chains with open boundary conditionsEl: 
The latter is obtained (in the AKLT limit) by taking 6 — > in our wavefunction. In the 
ladder the wavefunction is more general than just a RVB ansatz, but the effect of the 
quasifree boundary spins is the same. In this sense our results are related to the results 
of Hida0 who investigated the ladder with ferromagnetic coupling A on the rungs. This 
coupling connects two spins which in our antiferromagnetic ladder are not coupled directly 
but as nnn by two af bonds; from our calculations we see that this leads to the same effects. 

An alternative approach to arrive at a theoretical estimate for Jefj^ is to integrate out 
the spin degrees of freedom on the rungs between the boundary spins. For this purpose we 
write the general Hamiltonian with boundary spins 5*0 and >S'2p+i as 

H = Hiadderi^, 2...2p) + /Iq,! + /i2p,2p+l- (8) 

We eliminate /io,i and /i2p,2p+i to first order by a suitable canonical transformation and 
average over the eigenstates of the complete rungs forming the intermediate part of the 
ladder to end up with an effective low-energy Hamiltonian of the form 

H^^^ = Eo,±ip) + Jeff,±ip) So S^p+i. (9) 
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The explicit results for two and three intermediate rungs between the boundary spins are 

Eo^±{3) = -0.344, ^e/f+(3) = -0.2377, Jeff'^-i^) = 0.2248. (10) 

These results agree to the numerical data within 20%. For two isolated spins, higher 
order contributions to the canonical transformation will change the numbers towards the 
correct values but due to isotropy and to the fact that we are dealing with S = ^ the form of 
the Hamiltonian will not be affected. Although this approach gives only rough estimates, it 
serves to illustrate that for more than two unpaired spins higher order terms in the canonical 
transformation will lead to pair interactions between unpaired spins which are not nearest 
neighbors and to m-spin interactions (m > 4). The difference between the above analytical 
result and the numerical results should be considered as an indication that these higher spin 
interactions cannot be neglected at the outset and a description of the impure chain in terms 



of an effective pair Hamiltonian for unpaired spins as in ref. 12 needs justification. We will 
return to this question in section IV. 



Spin configurations in the presence of two unpaired spins 

We now discuss in more detail the spin configuration in the presence of two unpaired 
spins separated by a number of rungs larger than the correlation length. Then the spin 
projection will spread into the adjacent part of the ladder, i.e. into the rungs to the right 
and to the left of its site (or into the rungs on one side only in the case of edge spins) . In 
order to demonstrate this effect for an unpaired spin with a given projection (5*^) = ±| we 
consider the quantity {S^ _^), where m is the distance from the unpaired spin (m = denotes 
the unpaired spin) and ± distinguishes between spins on the same (+), resp. opposite (-), 
sublattice as the unpaired spin. When the two unpaired spins are sufficiently separated we 
obtain a clear picture of the situation in the neighborhood of a single unpaired spin coupling 
the two unpaired spins to S^ot = 1 and considering states with S^^^ = +1. Results from the 
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DMRG and Lanczos which illustrate this redistribuiton of the spin projection are shown in 
fig. ^; all data are of the following form: 

{^m,,a) = ^center 171 = 

= staii e'~ for |m| > 0, a = + 

I m| 

= ~ ^taii ^ ^ + for |m| > 0, a = - 



When the two unpaired spins are sufficiently far apart, the total spin projection +| has to 
be recovered, i.e. 

oo 2 

^center ~^ ^ ^ ~ 7^ 
m=l ^ 

for an unpaired spin in the bulk, and 

oo -j^ 

^center + "^P ~ 9 i^^) 

m=l ^ 

for an unpaired spin at the edge. (5p is a correction which approaches zero rapidly for p > 5. 
Quantitative results from the two approaches are: 
For an unpaired spin in the bulk (from DMRG): 

e = 3.1 s^,,,,, = 0.294 sL,z = 0.247, (12) 

and for an unpaired spin at the edge (from Lanczos, after corrections for the finite size of 
the system, i.e. the contributions of the second impurity): 

i = 3.1 ^Lnte. = 0.348 sl^^ = 0.308. (13) 

The amplitude s^^^i characterizes the redistribution of spin projection into the adjacent 
rungs and is seen to be rather large. The total excess spin on the positive sublattice is 0.894 
(bulk), resp. 1.117 (edge) and on the negative sublattice -0.792 (bulk), resp. -0.965 (edge). 
We notice that the asymptotic exponential behavior governed by the correlation length is 
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correct for all sites on the sublattice of the unpaired spin whereas it is approached within 
typically 5 rungs on the opposite sublattice. 

The numbers given above are consistent with the difference between the effective ex- 
change constants J^<^fl^ and J^''^ noticed above and give a microscopic understanding of this 
difference as is seen in the following way: For an unpaired spin in the bulk of the ladder 
the magnitudes of both s^^^^^ and s^^j^ are reduced by a factor of ~ .8 as compared to an 
unpaired spin at the ladder edge. Since the effective interaction of two unpaired spins will 
be determined by the tails of the spin distribution, Jg// will be smaller for unpaired spins 
in the bulk ladder by 

(z,bulk \ 2 
-0-64, (14) 

consistent with the value re^jc ~ -65 found from the analysis of the energy levels above. 

A redistribution of in the neighborhood of an unpaired edge spin is also obtained 
from the MP wave function above which gives 

Kenter ^ a' + ab ^ 0.422. (15) 

The behavior on the adjacent rungs is different in detail from the numerical results: 
Sp decays with alternating sign and purely exponential with different amplitudes on the 
leg of the unpaired edge spin (s^^-^ — + ab ^ 0.422) and on the opposite leg {s^^^i — 
—a/ {a + b) ^ —0.767). Again, the MP approach gives a reasonable qualitative picture but 
fails quantitatively. 

IV. MANY IMPURITIES 

It is natural to assume that the results obtained for two unpaired spins can be generalized 
to give results for the low energy properties of the general impure ladder with impurity 
concentration c by using pair exchange interactions between unpaired spins which are nearest 
neighbors only and determining sign and magnitude of these interactions from the results 



13 



in section III. For a random distribution of impurities this means that the impure ladder 
reduces to a spin chain with exchange constants which are random with respect to both sign 
and magnitude. This is the model used in ref. 0. 

In this section we will investigate this assumption more closely by comparing exact 
(DMRG) spectra for various configurations with four and six impurities to spectra obtained 
from the assumption above. We have restricted ourselves to an even number of impurities 
in order to facilitate the comparison between different impurity concentration (for an odd 
number of impurities the spin of the ground state will be half integer). A similar approach 
to the case of the spin-Peierls substance CuGeOs with impurities has been done in ref. |^. 

Before presenting these numerical results we want to discuss some aspects of the general 
impure ladder with antiferromagnetic exchange interactions only and L rungs, i.e. 2L sites, 
of which sites are occupied by nonmagnetic impurities. The ground state of the remaining 
2L — N spins S = j can take values Stot of total spin between and N/2, depending on 
the distribution of the impurities on the two sublattices. For a random distribution of the 
impurities the distribution of Stot is easily obtained from the theorem of Lieb and Mattis 
which states Stot = ^\^a — Nb\, where Na, Nb are the numbers of defects on the A and B 
sublattices respectively (for configurations with two defects on the same rung, which break 
the ladder, this gives Stot only for one out of a number of degenerate ground states). Thus 
for a random distribution of defects the probability g{Stot) for a ground state to have total 
spin Stot is 

for Stot = g{Stot) = (^^2) 

for 5„>0 ,(s„) = 2(^^^'^ J. (16) 

For ^ 1 we obtain from Stirlings formula the relative weight of states with Stot = and 
Stot = 1 (2^ is the total number of defect induced states) 



giStot = 0) 2 giStot = 1) 



(17) 



2^ V ttAT' 2^ V T^N' 

The ground state spin value with largest probability therefore is = 1. Of particular 
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interest are the average values of the ground state spin for a random distribution of ^ 1 



impurities. They are calculated using eq. (ITB) to give: 



(SL) = J (18) 

Thus the result obtained in ref. |T^ from a random walk argument is seen to be an exact con- 
sequence of the Lieb-Mattis theorem. For a discussion of the low temperature susceptibility 
in some given defect configuration we use 5*0 to denote the total spin of the ground state 
(with zero energy) and characterize the remaining states a = 1,2. ..am = 2^ — 1 by their 
total spin Sa and excitation energy Aq, > 0. The limiting susceptibility for low magnetic 
fields is then given by 

^ SkeT (25o + l) + Ea=i(25„ + l)e-^^«' ^ ^ 

with 

z{S^)=2S^{S^ + ^){S^ + l). (20) 

If the ground state of the impure ladder is not a singlet, i.e. 5*0 > 0, leading to z{So) > 0, a 
Curie susceptibility results in the low temperature limit. For So = 0, leading to z{So) = 0, 
the susceptibility shows an activated temperature dependence characteristic of a gapped 
system (Ai finite); if the low-lying states get dense with an asymptotic density of states 
p(e) ~ e~" the low temperature behavior of the susceptibility changes to x{T) ~ T~". In 
particular a temperature independent susceptibility results for a = 0, i.e. a constant density 
of states (this is the case of the 5" = | Heisenberg antiferromagnetic chain - actually it should 
be realized in this context that a concentration of c = ;| is sufficient to turn the ladder 
geometry into that of a single chain). All these deviations from a Curie behavior require a 
singlet ground state, which according to eq.(|T^ occurs with negligible weight in the limit 
of a macrosocpic system. We therefore conclude that the low temperature susceptibility of 
the impure ladder always follows a Curie law. This confirms from a different point of view 
the result a = 1 which has been obtained from a renormalization group approachEl. 
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When the hmiting value of the Curie constant for T — > is calculated on the basis of 
the above expressions, we find that due to the random positions of the defects the factor | 
(resulting from S'(S' + 1) for 5 = ^) is replaced by |. The experimentally interesting behavior 
at finite temperatures involves the transition between these two limiting cases, i.e. a change 
in Curie constant by a factor of 3. It will be determined to a large extent by the density of 
states. This is not easily accessible, and in the remainder of this section we present what 
can be learned from numerical calculations. 

In fig. ^ and table |T| we present a number of spectra with levels classified according to 
total spin 5* and corresponding results from the effective pair model. We have chosen the 
parameters under the following aspects: 

1. We want to control the effective model, 

2. we want to illustrate the filling of the gap with defect states, and 

3. we want to present the effects of the different sign combinations for defects on a given 
sequence of rungs. 

For the comparison between the complete ladder spectra as obtained from the DMRG 
and the spectra of the effective model we refer to table |T|. We see that the agreement is very 
satisfying, deviations are generally below 10% with a tendency of better accuracy for low 
energies. We note that all DMRG energy spectra overestimate the true energies. Further 
corrections will therefore rather improve the agreement between effective model and DMRG 
calculations. 

Since the spectra obtained by the two approaches agree so well, in fig. || DMRG spectra 
(full lines) only are given when these are available. However, the higher singlets and triplets 
were not accessible at reasonable computational expense as the respective states are already 
quite high-lying in the S^^^ = 0, 1 sectors. It is only in these cases that we have included 
into fig. 1^ the spectra obtained in the effective model (light gray lines). 
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In fig. I^(a-e) we show spectra for four unpaired spins with 4, 3 and 5 complete rungs 
in between. The possible combinations of the signs are seen to lead to quite different 
spectra and the general tendency to fill the regular ladder gap is evident. Figs. |^(a, e) 
result from each other qualitatively by an overall sign change. Although the magnitudes of 
the interactions are somewhat different (compare table ^ this is clearly evident in these 
spectra. We therefore present only three out of the remaining six spectra for this combination 
of distances in fig. ^(c - e). Fig. ^(f, g) presents two examples of symmetric configurations 
of four spins (here levels are characterized by the additional quantum number parity which 
we have not indicated). The six-defect spectrum shown in fig. |^(h) again is for a symmetric 
configuration. In particular these data illustrate, beyond related calculations, which have 
been done for the case of two impuritieJ^, how the density of states with S = 1 may increase 
at low energies with the number of impurities. However, a much larger number of impurities 
will be required to obtain a reliable numerical estimate for the low energy density of states 
Pie)- 



V. CONCLUSIONS 

We have demonstrated that interacting impurities drastically change the low-energy spec- 
trum of an antiferromagnetic spin ladder. A statistical distribution of impurities reduces 
the ladder, as far as its low energy spectrum is concerned (energy range of the pure lad- 
der gap), to a 5* = I chain with random interactions, a model which has recently been 
used to discuss the expected low temperature properties of the impure laddei§. From our 
calculations precise information on the parameters of the effective interactions is available. 
For large distances between defects the effective coupling strength is found to decay purely 
exponentially. We have found that it is sufficient to use two spin interactions to describe the 
spectrum of the effectively random chain to within 10%. The spectra which we have calcu- 
lated also illustrate the beginning of the process that a random distribution of impurities 
with concentration c produces a large number (of the order of 2'^^) of low- lying states which 
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will fill the ladder gap as observed experimentallyO. Applying the Lieb-Mattis theorem to 
the impure ladder with antiferromagnetic interactions only we argue that for T ^ the 
susceptibility is characterized by a 1/T divergence with finite Curie constant. Our numer- 
ical results, however, are not sufficiently accurate to draw conclusions about the change in 
Curie constant with increasing temperature. 

From our results it becomes also clear that interactions between unpaired spins which are 
not nearest neighbors and m— spin interactions are present and cannot be neglected in the 
effective model at the outset. They have, however, turned out not quantitatively relevant 
for the low-energy spectra which we have computed. 

The effects of impurities on the spin distribution which we have presented show some 
similarity to what has been discussed in the context of impure Haldane chain^ and it is useful 
to qualitatively compare these two phenomena once more: In both systems the unpaired 
spin at the position of the impurity shares its magnetic moment with the adjacent part of the 
system - and the degree of this mixing, i.e. its spatial extent, is determined by the correlation 
length. In the Haldane case, one deals with a physically different magnetic impurity spin; it 
interacts with that part of the unpaired boundary spin which remains localized, leading to 
the experimentally observed splitting of the ESR spectra. In the spin ladder a nonmagnetic 
impurity is present and the corresponding unpaired spin shares its magnetic moment with 
the adjacent parts of the ladder; the tails of this magnetic moment introduce an effective 
interaction leading to the experimentally observable splitting of the low-energy spectrum. 
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FIGURES 

FIG. 1. Structure of a spin ladder with three impurities leading to unpaired spins in cis resp. 
trans configuration. 

FIG. 2. Effective exchange interactions for two unpaired spins at the edges and in the bulk of 
a spin ladder. Exchange energies for effective ferromagnetic and antiferromagnetic interactions are 
plotted separately. 

FIG. 3. Low-lying energy levels for ladders with periodic boundary conditions and two unpaired 
spins with a distance of pi resp. p2 complete rungs. Spectra for unpaired spins with effective 
antiferromagnetic, resp. ferromagnetic, interactions are shown separately, (a, b) pi + P2 = 12 
rungs, (c, d) pi + p2 = 11. 

FIG. 4. Distribution of the excess z-component of spin (magnitude 1/2) resulting from the 
unpaired spin opposite to an impurity into the adjacent rungs of the ladder (the positive and 
negative contributions for each rung p are alternating between the two legs). 

FIG. 5. Low energy spectra of the ladder with four and six impurities from DMRG (black 
lines) and from the effective model (light gray lines, only when DMRG data are not available) 
for different configurations characterized by the sign (F, AF) and the number of complete rungs 
between the unpaired spins as indicated in the graphs, (a) - (g) four unpaired spins: (a) F4 - F3 
- F5, (b) AF4 - F3 - AF5, (c) F4 - AF3 - AF5, (d) F4 - F3 - AF5, (e) AF4 - AF3 - AF5, (f) F5 - 
AF4 - F5, (g) F3 - AF4 - F3; (h) six unpaired spins: AF5 - F4 - AF6 - F4 - AF5. 
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TABLES 

TABLE L Energies calculated by the DMRG before (first line) and after one (second line) and 
two (third line) applications of the finite size algorithm. 

TABLE IL Effective exchange energies for two unpaired spins at the edge and in the bulk of 
a ladder with p rungs in between. Two energies are given for each p, corresponding to the two 
positions on each rung. 

TABLE in. Low energy spectra of the ladder with four and six impurities for different config- 
urations as in fig.|5|: comparison of DMRG and effective Hamiltonian results. 



21 



Fig.1 



Trans configuration 



Cis configuration 



Fig. 2 



G--OAFM (boundary) 
G- "OAFM (bulk) 
A--AFM (boundary) 
A - -A^'^ (bulk) 




8 



10 



P 



Table 2 



Jeff (ED) 


rungs (n) 


antiferromagnetic 


ferromagnetic 


2 


0.4134342076 


0.3332389325 


3 


0.2808007550 


0.2790132929 


4 


0.1944368705 


0.2010576140 


5 


0.1379125645 


0.1459181428 


6 


0.0990135730 


0.1051339920 


7 


0.0717149472 


0.0759464504 


8 


0.0522047020 


0.0549482466 


9 


0.0381181864 


0.0398361954 


10 


0.0278774960 


0.0289263601 


11 


0.0204037286 


0.0210324667 


12 


0.0149374382 


0.0153089483 


Jeff (DMRG) 


rungs (n) 


antiferromagnetic 


ferromagnetic 





0.471 




1 


0.398 


0.249 


2 


0.255 


0.170 


3 


0.181 


0.154 


4 


0.126 


0.115 


5 


0.090 


0.087 


6 


0.065 


0.064 


7 


0.047 


0.047 


8 


0.034 


0.035 


9 


0.025 


0.025 


10 


0.018 





Fig. 3.a 



-13.0 - 



c 

0) 



-13.5 - 



-14.0 - 



-14.5 




Fig. 3.b 



0) 

c 

0) 



-12.8 - 



-13.1 - 



-13.2 - 



-13.4 - 



-13.7 - 



-13.8 - 



-14.1 - 



-14.2 




Fig. 3.C 



-12.0 - 



O) 

c 

0) 



-12.4 - 



-12.8 - 



-13.2 





1 ■ 1 


1 


1 


1 ■ 1 ■ 

Osinglet 










Atriplet 




A 


A 


O 


A 






A 


A 


A 




A 










A 


A 


A 


A 






O 


O 






O 

1 . 1 


1 


1 





3 4 
Pi 



Fig. 3.d 



-12.0 - 



O) 

c 

0) 



-12.4 - 



-12.8 - 



-13.2 



■ 1 


1 


1 


1 


1 ■ 1 ■ 
Osinglet 










Atriplet 








o 


- 

o 


A 




A 


A 


O 










O 












O 


o 






A 


A 


A 


Q 
A 




) 1 


2 


3 


4 


5 6 7 



P1 



Table 3 



A(4) F(3) A(5) 


A(4) A(3) A(5) 


gtot 


ladder 


eff. model 


gtot 


ladder 


eff. model 





0.000 


0.000 





0.000 


0.000 


1 


0.079 


0.07776 


1 


0.047 


0.04730 


1 


0.131 


0.12681 


1 


0.197 


0.19146 


2 


0.198 


0.19196 





0.204 


0.20080 


1 




0.30932 


1 


0.269 


0.26093 







0.32193 


2 


0.310 


0.29890 


F(4) A(3) A(5) 


F(4) F(3) A(5) 


1 


0.000 


0.000 


1 


0.000 


0.000 





0.030 


0.02582 


2 


0.075 


0.07296 


1 


0.156 


0.15188 





0.094 


0.09089 


2 


0.208 


0.20117 


1 


0.134 


0.12860 







0.22053 







0.23404 


1 


0.302 


0.29564 


1 




0.26930 


F(4) F(3) F(5) 


F(5) A(4) F(5) 


2 


0.000 


0.000 





0.000 


0.000 


1 


0.036 


0.03474 


1 


0.010 


0.010 





0.094 


0.08660 


2 


0.109 


0.107 


1 


0.108 


0.10041 


1 


0.168 


0.164 


1 




0.22085 


1 


0.201 


0.194 







0.26940 





0.273 


0.261 


F(3) A(4) F(3) 


A(5) F(4) A(6) F(4) A(5) 





0.000 


0.000 





0.000 


0.000 


1 


0.017 


0.017 


1 


0.051 


0.049 


2 


0.109 


0.104 


1 


0.088 


0.081 


1 


0.225 


0.216 


1 


0.095 


0.091 


1 


0.271 


0.257 


2 


0.128 


0.126 





0.403 


0.387 


2 


0.148 


0.133 













0.161 








1 




0.166 








2 


0.196 


0.182 








3 


0.216 


0.202 








1 




0.219 








1 




0.223 













0.226 













0.235 








2 


0.298 


0.281 








1 




0.309 








2 


0.325 


0.311 













0.398 








1 




0.401 



0.40 



0.20 - 



0.00 



Fig. 4 



-0.20 

-7 -5 -3 -1 1 



Fig. 5 



C 
0) 



0.4 - 



0.2 - 



0.0 - 



S=0 S=2 
S=1 



c 

0) 



0.4 - 



0.2 - — _ 



0.0 I- — 

S=0 S=2 
S=1 



S=0 S=2 
S=1 



S=0 S=2 
S=1 



S=0 S=2 
S=1 



S=0 S=2 
S=1 



S=0 S=2 
S=1 



S=0 S=2 

S=1 S=3 



Table 1 



S' = 


S' = 1 


M = 40 


M = 100 


M = 40 


M = 100 


-65.495857 


-65.505486 


-65.089740 


-65.103673 


-65.518368 


-65.518391 


-65.118436 


-65.120279 


-65.518388 


-65.518391 


-65.120622 


-65.120652 



